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We provide a comprehensive model-independent analysis of rare decays 
involving the b — > s transition to put constraints on dimension-six AF = 1 
effective operators. The constraints are derived from all the available up-to- 
date experimental data from the B-factories, CDF and LHCb. The implica- 
tions and future prospects for observables in b — > si + i~ and b — > svv transi- 
tions in view of improved measurements are also investigated. The present 
work updates and generalises previous studies providing, at the same time, 
a useful tool to test the flavour structure of any theory beyond the SM. 



1. Introduction 



The CKM description of flavour and CP violation in the Standard Model (SM) has 
been extremely successful in describing the data from high-precision flavour experiments, 
among them the B factories and the Tevatron, and is also in agreement with recent 
measurements of rare B decays at LHC. This is a remarkable fact, given that the origin 
of flavour, i.e. the origin of the replication of fermion species and of the hierarchies in 
their masses and mixing is still a complete mystery. If one supposes new physics (NP) to 
be at work not too far above the electroweak scale, as is implied by the gauge hierarchy 
problem, this fact is even more surprising, since NP coupling to the SM generally leads 
to modifications of flavour violation, in particular in flavour-changing neutral current 
(FCNC) processes, which are suppressed by the GIM mechanism and arise only at the 
loop level in the SM. Even under the most restrictive of assumptions that can be made on 
the flavour sector of a NP theory, namely assuming all flavour violation to be governed 
by the SM Yukawa couplings (the hypothesis of Minimal Flavour Violation (TJ[2j , MFV) , 
many models predict significant deviations in FCNC observables. Moreover, many of the 
flavour-violating couplings are rather poorly constrained as yet and still allow for sizable 
NP contributions. In AB = 1 transitions, this is particularly true for NP contributions 
with a different CP phase or chirality with respect to the SM contribution. In that case, 
large deviations in observables not measured yet are still easily possible. Quantifying 
this statement is one of the main goals of this paper. 

Experimentally, the prospects to improve constraints on flavour violating couplings 
are excellent: Today, the LHCb experiment has a high sensitivity to exclusive hadronic, 
semi-leptonic and leptonic B and B s decays |3|J^] In the mid-term future, two next- 
generation B factories will allow also the measurement of inclusive rare decays and 
decays with neutrinos in the final state (4j[5]. 

Since FCNC processes test NP indirectly, through quantum corrections induced by 
heavy particles, the impact of NP on low-energy observables like branching ratios or 
asymmetries can be summarised by their modification of Wilson coefficients of local, 
non-renormalizable operators. These short-distance coefficients can be constrained on 
a completely model-independent basis by measuring FCNC observables and these con- 
straints can in turn be used to constrain individual NP models. 

In this paper, we concentrate on AB = AS = 1 processes, i.e. rare decays with 
a b — > s transition. We use up-to-date experimental constraints - in particular, we 
include the recent measurement at LHCb of angular observables in B — > K* fi + yT [6j, 
which is the most precise to date - to put model-independent constraints on the Wilson 
coefficients Similar constraints have been considered in the literature before, e.g. in the 
context of MFV 19 - 11 , magnetic penguin operators 12 , SM operators 13 -15 or generic 



NP 16-18 . We generalise these studies by considering SM operators, their chirality- 



1 In the case of leptonic decays, also ATLAS and CMS are competitive. 

2 Instead, we do not consider here observables related to the decay B — > K£ + £~ , since their current 
experimental resolutions are rather poor. Including them would affect our results only in a negligible 
way. However, once LHCb data on B — ¥ K£ + £~ observables will become available, it will be important 
to include them given their high NP sensitivity [7l[8l- 
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nipped counterparts and generic CP violation and considering also the case where they 
are all simultaneously present instead of considering only a pair at a time. After imposing 
the above constraints, we investigate the room left for NP in observables which have not 
been measured yet, like the branching ratios of B s — > and B s — > t + t~, CP 

asymmetries in B — > K*fj, + fj,~ and observables in b — > suT> decays. 

The paper is organised as follows. In section [2| we specify the effective Hamiltonian 
for b — >• s transitions and discuss all the observables relevant for our study. In section [3j 
we present our numerical results for the model-independent constraints on the Wilson 
coefficients. In section [4j we consider the constraints in the more restrictive case of semi- 
leptonic operators generated dominantly by modified Z couplings, which is the case in 
many models beyond the SM. This allows us to correlate b — > s£ + £~ processes to b — > svv 
processes. Our main findings are summarised in section [5} 

2. Observables in b — > S7 and b — >• si + t~ decays 

The effective Hamiltonian relevant for b — > s'j and b — > s£ + £~ transitions is given 
by [19|[20 

n eS = -^ V M^ + + h - c - • C 1 ) 

i 

The operators 0{ that are most sensitive to NP effects are 

7 = ^(Kt^Pr^F^, Os = 9 -^{sa^T a P R b)G^\ 

9 = (s^P L b) (If I) , Ova = (s^P L b) {hV) , 

O s = m b (sP R b){££), Op = m b (sP R b)(£ l5 £) , (2) 

where nib denotes the running b quark mass in the MS scheme and Pl,r = (1 -F 7s)/2- 
The corresponding operators 0\ are obtained from the operators Oi via the replacement 
Pl P R . Numerical values of the Wilson coefficients and the relation of the coefficients 
at the matching scale to the effective low-energy coefficients are discussed in appendix [Bj 
Note that we assume and cjg to be independent of the lepton flavour, but Cg 

and Cp^ to be proportional to the lepton Yukawa couplings (this assumption will become 
relevant when we discuss the relation between B s — > and B s — > r + r _ ). Moreover, 

we also assume lepton flavour conservation, which is an excellent approximation for our 
purposes, given the stringent experimental bounds on lepton flavour violating processes. 

We will now discuss the observables in processes sensitive to this effective Hamiltonian, 
which can be used to constrain new physics. 

2.1. B s -> £+£- 

In the SM, the B s — > fi + fi~ decay is strongly helicity suppressed and among the rarest 
FCNC decays. Using directly the value of the B s meson decay constant from the lattice, 
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f Bs = (250 ± 12) MeV |2l| we get the following SM prediction^ 



BR(B S -»• fi + fi~) S M = (3.7 ± 0.4) x 10~ 9 . (3) 

The LHCb and CMS collaborations have set a combined upper bound on the branching 
ratio of [24fj26 

BR(B S -> m + A*~)lhc < 1-1 x 10- 8 (4) 

at 95% confidence leveljf] 

In a generic NP model, the branching ratio of BR(i? s — > is given by 



BR(B S -> / u+^-) SM y m B 

where 



l-sf i-z^ KIT. (5) 



(Cs - C 's) P _ (Cp - Op) ^ (Cio - Ci ) 



" 2m, |C? M | ' 2m M Cf M Cf M ' 1 j 

The important feature of B s — > /i + as a probe of NP is that it is among the very 
few b — > s decays that are strongly sensitive to scalar and pseudoscalar operators. In 
models where such operators are sizable, the branching ratio can easily saturate the 
experimental limit. In section [3j we will use measurements of other b — > s processes 
to constrain the Wilson coefficients c{q, which then allows us to predict the maximum 
allowed size of BR(£? S — > /i + /i~) in the absence of (pseudo) scalar currents. 

The process B s — > t + t~ , is governed by Wilson coefficients analogous to B s — > 
and its branching ratio is given by eqs. ^ and Q with the appropriate replacement 

fjb — ^ T. 

2.2. b ->■ sj 

The experimental data on the branching ratio of the inclusive B — > X s ^ decay |28| 

BR(B -> X s7 ) cxp = (3.55 ± 0.26) x 10 -4 , (7) 

and the corresponding NNLO SM prediction [29[|30| 

BR(B -»• X s7 ) S m = (3.15 ± 0.23) x 10" 4 , (8) 

show good agreement. As the BR(i? —> X s j) is highly sensitive to NP contributions 
to the Wilson coefficients Cj , this agreement leads to severe constraints on the flavour 



3 Using the remarkably precise value for fs s obtained very recently in [22]: /b s = (225 ± 4) MeV, 
we obtain BR(_B S — > /i + /i~)sM = (3.0 ± 0.2) x 10 -9 . This result is as precise as the value that is 
obtained by assuming AM S free of NP and using its me asur ement to reduce the theory uncertainty in 
B s — > n + fi~ arising from the B s meson decay constant [23]. However, to be conservative, we use ([3]) 
in our analysis. 

4 We mention that the CDF collaboration found an excess of candidates for B s — > decays, which 



have been used to determine 27 BR(_B S — > /i + fi )cdf = (l.8_o!g) x 10 
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sectors of many NP models. In our numerical analysis, we use the expression for the 



branching ratio reported in 31 , rescaled to the SM prediction of |29|, and assume the 



uncertainty in the SM prediction as relative error on the theory prediction. 

Another interesting observable that probes the b — > sj transition is the time-dependent 
CP asymmetry in the exclusive B& — > K*(—> Kg7r°)j decay 



32 



35 



T{B°(t) -»■ if* u 7 ) - T(B°{t) if* u 7 ) 
T(B°(t) -> K*°j) + T(B°(t) K*°-f) 



Sk**, sm(AM d t) - C K *-y cos{AM d t) . (9) 



The coefficient Sk*^ is highly sensitive to right handed currents as at leading order it 
vanishes for C 7 — > 0. As a consequence, SM contributions to Sk*j are suppressed by 



iris/mi, or Aqcd/^J) [33], resulting in a very small SM prediction 35 



c-SM 



(-2.3 ± 1.6)% 



(10) 



Experimental evidence for a large Sk*-v would be a clear indication of NP effects through 



right handed currents. On the experimental side one has presently 28,36,37 



QCXp 



-0.16 ±0.22 



fill 



and the prospects to improve this measurement significantly at next generation B facto- 
ries are excellent [38]. In our numerical analysis we use the LO expression for Sk*~j [34] 



Sk* 



\Cr\ 2 + \C' 7 \ 



,-Im ( e 



-i<t>d 



c 7 c 7 



(12) 



that leads to accurate predictions in presence of NP. In the above expression for Sk*^, 
sin(^) = S^Ks i s t ne phase of the B d mixing amplitude and the Wilson coefficients 
are evaluated at the scale \x = m&. Using directly the experimental value for S^k s = 
0.67 ± 0.02 [28], we automatically capture possible NP effects in B d mixing. In our NP 



analysis, we assume the theory uncertainty to be equal to the SM uncertainty in (10). 

Another observable that is in principle sensitive to CP violating effects in the 6—7-57 
transition is Acp(b —> 57), the direct CP asymmetry in the B — > Xg'j decay 39,40 . In 



contrast to the observables discussed so far, it is also highly sensitive to NP contributions 



to the Wilson coefficients Cg'^ . However, as shown in 
sj) is dominated by long-distance contributions and 
difficult to predict this observable reliably in the context of NP scenarios. We therefore 
do not consider Acp(b — > S7) in our analysis. We also do not consider the isospin 
asymmetry in B — > if* 7, as large hadronic uncertainties strongly limit the constraining 
power of this observable. 
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, the SM prediction for Ac?(b 



arge hadronic uncertainties make it 



2.3. B -)■ X s £ + i~ 

We consider the inclusive B — > X s £ + l~ decay in two different regions of the dilepton 
invariant mass. The low q 2 region with 1 GeV 2 < q 2 < 6 GeV 2 and the high q 2 region 
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Table 1: Dictionary between different notations for the -B — )• K*fi + [i~ observables and 
Wilson coefficients they are most sensitive to (the sensitivity to ci is only 
present at low q 2 ). 



with q 2 > 14.4 GeV 2 . Averaging the available results from BaBar |42| and Belle 43 one 
finds the following averages for the branching ratios in the two regions 

BR(B -> X s £ + r)™ p 6] = (1.63 ±0.50) 10" 6 , BR(B -> X s l + rf^ Ai = (4.3 ±1.2) 10" 7 . 

(13) 



These results should be compared to the SM predictions 119, 44-51 



BR(B 



(1.59 ±0.11) 10~ b , BR(B ->• X. 



)>14:A 



(2.3 ±0.7) 10" 7 . 

(14) 

While the low q 2 values are in perfect agreement, the SM prediction in the high q 2 region 
is on the low side of the experimental result. We remark that the theory prediction in 
the low q 2 region quoted above does not include effects from the experimental cut on the 
hadronic final state. Such effects can be as large as 10% 52,53 . Enlarging the theory 



error correspondingly however would effect our results only in a minor way, given the 
huge experimental uncertainty in BR(S X s l + £~). 

The branching ratio in the high q 2 region is mainly sensitive to NP contributions to 



the Wilson coefficients and C)q , while the branching ratio in the low q 2 region also 
depends strongly on Cj\ In our numerical analysis, we use the expressions given in |ll| , 
adjusting them to take into account also the primed Wilson coefficients, and treat the 
uncertainties in the SM predictions as relative errors on the theory predictions. 

In principle, another interesting observable to constrain NP would be the forward- 
backward asymmetry in B — > X s £ + £~ . However, since it has not been measured yet, we 
do not include it in our analysis. A fully inclusive measurement is probably not feasible 
at LHCb and it will only be possible at next generation B factories [5,38,54 . 



2.4. B ->■ K 



The angular distribution of the exclusive B — > K* (— > K~ tt + ) /x" 1 " \x~ decay gives access to 
many observables potentially sensitive to NP [14, 15,17,31] 55, 56,61-65]. By means of its 
charge conjugated mode B — > K*°(— > K + ir~)[i + pT , which can be distinguished from the 
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former simply by the meson charges, this decay allows a straightforward measurement 
of CP asymmetries. 

Neglecting scalar operator contributions (which are strongly constrained by B s — > 
and lepton mass effects (which is a very good approximation for electrons and 
muons even if effects from collinear QED logarithms are taken into account |48.|), the 
full set of observables accessible in the angular distribution of the decay and its CP- 
conjugate is given by 9+9 angular coefficients Ii(q 2 ) and h{q 2 ), which are functions of 
the dilepton invariant mass q 2 . While the overall normalization of the angular coefficients 
is subject to considerable uncertainties, theoretically cleaner observables are obtained by 
normalizing them to the total invariant mass distribution. Furthermore, it makes sense 
to separate the observables into CP asymmetries Ai and CP-averaged ones Si. One thus 
arrives at 



55 



4=H + f.)/^, A -=^- I ')/ dJ W 1 - W 

We will also consider observables integrated in a q 2 range, defined as 

<Si>m = (/V (/* + /*)) / (/V^J^) , (16) 

and analogously for {Ai}. 

For all observables one has to distinguish, both theoretically and experimentally, be- 
tween the kinematical region where the dilepton invariant mass is below the charmonium 
resonances (low q 2 or large recoil region) and the region above (high q 2 or low recoil re- 
gion). The intermediate region is of no interest to probe NP, as the cc resonances 
dominate the short distance rate by two orders of magnitude. 

At low q 2 , the observables are sensitive to all the Wilson coefficients \ 10 . Among 
the CP asymmetries^J the most promising ones are then the T-odd CP asymmetries 
Ay, A§ and Ag, which are not suppressed by small strong phases [17] ■ At high q 2 , the 
contributions of the magnetic penguin operators Cy are suppressed, which in turn allows 
a cleaner sensitivity to the semi-leptonic operators. The CP asymmetries reduce to three 
independent ones, which are however T-even and therefore suppressed by small strong 



phases even beyond the SM 15 . Among the CP-averaged angular coefficients, two have 



already been measured (6 57 - 60 : the forward-backward asymmetry Ap~B and the K* 
longitudinal polarisation fraction Fj_,. Recently, the CDF collaboration also published 
first bounds on S3 and Ag [60| . A promising observable in the early phase of LHC is the 
observable S5 [67] . 

Since different notations and conventions exist for the numerous B — > K*£ + l~ observ- 
ables, in table [T] we provide a dictionary between the notation used in this work and a 
selection of other theory and experimental papers. It also lists the Wilson coefficients 



5 We do not consider the CP asymmetries A^^ s and A$ defined in 66 since the former is suppressed 
by a small strong phase even beyond the SM and the latter is normalized to the quantity Is + Is, 
which is zero at LO even beyond the SM and afflicted with considerable uncertainty. 
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which, if modified by NP, would have the biggest impact on the observable in question. 
In the case of C7 and CL this sensitivity is only present at low q 2 . 

The main challenge in the theoretical prediction of the B — > K*l + £~ observables is 
given on the one hand by the B — > K* form factors; on the other by non-factorisable 
effects^] At low q 2 , QCD factorisation can be used in the heavy quark limit, which 
reduces the number of independent form factors from 7 to 2 and allows a systematic cal- 



culation of non-factorizable corrections [69 70 . The remaining theoretical uncertainties 



then reside in phenomenological parameters like meson distribution amplitudes, in the 
form factors themselves, as well as in possible corrections of higher order in the ratio 
^■QCT>/ m b- Instead of using the two form factors in the heavy quark limit, we use the full 
set of seven form factors calculated by QCD sum rules on the light cone (LCSR), using 
the results of [55,71). This approach has two advantages. First, using the full set of form 



factors takes into account an important source of power suppressed corrections at low 
q 2 . Second, the correlated uncertainties between the different form factors obtained from 
the LCSR calculation leads to a strongly reduced form factor uncertainty on observables 
involving ratios of form factors. 

At high q 2 , QCD factorization and LCSR methods are not applicable. For the form 
factors, lacking predictions from lattice QCD, one currently has to rely on extrapolations 
of low-q 2 calculations, which introduce considerable uncertainty. For the estimation of 
non-factorizable corrections, an operator product expansion in powers of l/yip can be 



used 72 73 and in Ref. [73] it has been argued that non-perturbative corrections not 
accounted for by the form factors are of the order of only a few percent. We do take 
into account non-factorizable corrections proportional to form factors at 0(a s ) both at 
low and high q 2 [45)[5l|[69)[70)[74] 

For our numerical analysis, a description of our treatment of theory uncertainties in 
the B — > K*n + [.i~ observables is in order. For both high and low q 2 , we take into account 
parametric uncertainties, varying the ratio m c /mb from 0.25 to 0.33, the renormalization 



scale from 4.0 to 5.6 GeV 55 and the CKM angle 7 by ±11° 75 . At low q 2 , as mentioned 
above, we make use of the LCSR calculation of all 7 form factors and vary the LCSR 
parameters as discussed in Ref. |55| . To be conservative, we add an additional real scale 
factor with an uncertainty of 10% to each of the transversity amplitudes to account 
for possible additional power suppressed corrections. The branching ratio is the only 
observable that is sensitive to the overall normalization of the form factors. Since LCSR 
only give predictions for the B meson decay constant times a form factor, we add an 



additional relative uncertainty of twice the uncertainty of fs = (205 ± 12) MeV 21 
to the branching ratio. Unlike in [55], we do not use the data on B — > K*j to fix 
the form factor normalization, since we allow for NP also in B — > K*j. We add all 
the individual uncertainties in quadrature. At high q 2 , we use the extrapolated form 
factors of Ref. [76] . In the Simplified Series Expansion used there, each form factor 
depends on two parameters fitted to the low-g 2 LCSR calculation. We estimate the 
form factor uncertainty by varying all 14 fit parameters separately, i.e. considering the 
uncertainties of the individual form factors as uncorrelated to each other, and add the 



6 For a recent discussion of uncertainties in the low q 2 region see also 68 
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resulting errors in quadrature. We consider this approach to be conservative. In view of 
the resulting sizable form factor uncertainties, at high q 2 we do not consider additional 
uncertainties due to power corrections or duality violation, which should amount to only 
a few percent [73] and are therefore numerically irrelevant. 

Figure [T] shows the predictions for Fl and ^4fb with our error estimates at low and 



high q and compares them to the experimental data from Belle 58 , CDF |77| and 



LHCb |6j. We do not show the data from BaBar [57], since they are given in large bins 
that include q 2 regions which are under poor theoretical control. 



3. Model-independent constraints on Wilson coefficients 

In a vast class of models beyond the SM, all NP effects in the observables listed in table [3] 
and discussed in section [2] are described by a modification of the Wilson coefficients 
Cj g g 10 s P at a matching scale, typically of the order of heavy particles contributing to 
the FCNC processes. For definiteness, we will consider in the following constraints on 
Wilson coefficients evaluated at the scale fih = 160 GeV. The values at any other scale 
below the matching scale can be obtained straightforwardly using the renormalization 



group 1 78 , 79 , see also appendix B 



Up to subleading contributions, the coefficients enter the processes of interest 
only via renormalization group running by the mixing of the operators and 0g ■ 
Therefore, we will not consider constraints on Cg in the following and keep in mind 
that, in the presence of Cg , the constraints on C^ NP we present can be understood as 
constraints on the combination (C^ NP + 0.16Cg'^ NP ) at [a^ (see appendix |b| . Among 
the observables we consider, the scalar and pseudoscalar coefficients C^ p can only affect 
the branching ratios of B s — > [i + and B s — > t + t~ in a significant way. Thus, we will 
disregard also these coefficients and instead use the constraints on the remaining Wilson 
coefficients to give a prediction for the maximum possible sizes of BR(i? s — > and 
BR(i? s —7- t + t~) in the absence of (pseudo) scalar operators. We are thus left with the 
6, potentially complex, Wilson coefficients \ 10 . 

To obtain constraints on the Wilson coefficients, we construct a x 2 function, which is a 
function of Wilson coefficients C and contains the theory predictions for the observables 
Of* and the experimental central values 0^ xp as well as the corresponding uncertainties 
(which we assume to be Gaussian), 

2 . _ (or-o?(c)Y 

y (C) = > — — (17) 

YK cxp ) 2 + (^)) 2 ' 

We write the theory uncertainty as a function of Wilson coefficients since, as discussed 
in section [2] it is a relative error for some observables and for some others, such as 
the B — > K* n + [i~ angular coefficients, can even be a non-trivial function of Wilson 
coefficients. 
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gure 1: Theory prediction (blue bands), binned theory prediction (purple horizon- 
tal bands) and experimental measurements from LHCb (red dots, thick er- 
ror bars), CDF (gray squares, dashed error bars) and Belle (gray triangles, 
solid error bars) for ^4fb and Fl in B — > K* ji + fi~ for the two theoretically 
controllable q 2 regions. 
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10 6 xBR(B X s £ + £-) [lj6] 


1.63 ±0.50 
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10 7 xBR(5 if*^+£-) [1)6] 


1.71 ±0.22 


6 
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2.28 ±0.63 
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1.11 ±0.13 
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77 
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10 7 xBR(B -> K*^ + r) [16il9] 


1.35 ±0.15 
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77 




1.34 ±0.51 
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0.28 ±0.09 
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(-0.3 ±1.1) 10~ 2 




(A 9 )(s^^+r) [li6] 
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60 
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Table 2: Experimental averages and SM predictions for the observables used in the fit. 



In table [3] we summarise the SM predictions as well as the experimental values of 
the observables that we use in the x 2 function. To obtain the experimental values we 
perform weighted averages of the available measurements, symmetrising the errors using 
the prescription of ref. |80| 



3.1. Impact of observables on pairs of Wilson coefficients 

Using the \ 2 function, we can obtain constraints on the real or imaginary parts of a pair of 
Wilson coefficients, or in the complex plane of a single Wilson coefficient. This approach, 
which is reminiscent of the CKM constraints in the p-fj plane, has the advantage that 
it allows to transparently show the impact of individual observables on the constraints. 
A similar approach has been used e.g. in 14,15,17] for the SM Wilson coefficients and 



m 



12 for Cj vs. C 7 . Moreover, in some cases a pair of Wilson coefficients (or a single 



complex coefficient) captures already the dominant NP effect in certain scenarios, in 
which case such plots become particularly useful. For example, in the MSSM with MFV 
and flavour blind phases |81| , in effective SUSY with flavour blind phases |82| and in 
effective SUSY with a U(2) symmetry 83,84], NP effects in AB = AS = 1 processes 
arise almost exclusively through complex contributions to C-j (and Cg)- In MFV models 
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with dominance of Z penguins and without new sources of CP violation, only the real 
parts of C7 and C\q are relevant as will be discussed in section |4j 

The resulting plots are shown in figure [2j In these plots, the dark and light red regions 
show the 1 and 2a best fit regions (contours of xiot ~ Xtot mm = 1 or 4) , while the shaded 
regions in different colours show the 2a allowed regions (contours of x 2 ~ X m in = 4) from 
B — > K*^L + fi~ at low q 2 (blue), B ->■ K* fi + yT at high q 2 (green), B ->■ X s £ + £~ (brown), 
BR(i? — > X s ^y) (yellow) and B —¥ if* 7 (purple). We only show the observables that give 
relevant constraints. 

We make several observations. 

• At the 95% C.L., all best fit regions are compatible with the SM. 

• In the complex C-j plane, which is relevant e.g. for the models with flavour blind 
phases mentioned above, the inclusive and exclusive b —> s£ + £~ observables - in 
particular the measurement of BR(S — > X s £ + £~) at low q 2 13 and the LHCb 



measurement of ^4fb a t low q - exclude a sign-flip in the low-energy Cf that 
would be allowed by BR(S — > X s ~f) and that was favoured by Belle data on 



^4fb 58 . We observe that an imaginary part of CV as large as |Im(CV)| < 0.7 is 



still allowed in this scenario. 

In the presence of Cf, the current data on the time-dependent CP asymmetry in 
B — > K*^ gives already an important constraint. 

Similar to the complex C7 plane, in the Re(C7)-Re(C 7 ) plane a sign- flip in Cf is 
excluded by the inclusive and exclusive b — > s£ + £~ observables. Imposing also the 
constraint from Sk*^ leaves only two disjoint regions at 95% C.L.: one around the 
SM point and a second, less favoured one with a large Re(C 7 ) ~ 0.5. 

In the complex C9 and C10 planes, sign flips of the real parts are excluded but 
imaginary parts as large as |Im(Cg)|, |Im(Cio)| ^5 3 are still allowed. 

The slight preference towards non-SM values of the Wilson coefficients in the com- 
plex Cg and C'i planes results from the tension between SM prediction and ex- 
perimental data on BR(i? — > X s £ + £~) in the high q 2 region as well as the tension 
between SM and experimental data on Fl in B — > K*fi + [i~ at low q 2 . The tension 
in Fl is also the reason for the slight preference for a non-SM value in the complex 
C' 7 plane. 

High-q 2 data on6-> K*£ + £~ are competitive with and coplementary to the low- 
q 2 ones. In particular, they are crucial to exclude a sign flip in Re(Cio) as well as 
simultaneous sign flips in C10 and C[ or C7 and Cg. 

In the Re(Cg)-Re(Cio) plane, a simultaneous sign flip in Cg and C10 would be 
allowed by high q 2 data on B — > K*£ + £~ but is excluded by the new Ap^ mea- 
surement at low q 2 from LHCb. 
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Figure 2: Individual 2a constraints on pairs of Wilson coefficients from B — > K* n + \i~ 
at low q 2 (blue) and high q 2 (green), B -> X s £ + £~ (brown), BR(B -)■ X s j) 
(yellow), B — > K*j (purple) and combined 1 and 2a constraints (red). 
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We stress that the above conclusions hold if one allows only the two quantities shown to 
be non-zero. In many NP models, several of the Wilson coefficients will deviate from the 
SM, which renders some or all of these constraints ineffective. The more general case 
calls for a global fit of all Wilson coefficients. 

3.2. Global fit of Wilson coefficients 

While the constraints discussed above are useful to display the constraining power of 
individual observables, they are not suited to put constraints on Wilson coefficients in 
models where more than two real or more than one complex coefficient is relevant, since 
cancellations can easily occur that render some of the constraints ineffective. Therefore, 
we now present constraints on Wilson coefficients varying not only 2 but all (or a sub- 
set) of the 6 complex coefficients Cj \ 10 . To cope with the large dimensionality of the 
parameter space, we perform a Markov Chain Monte Carlo (MCMC) analysis. Details 
on the statistical approach are given in appendix [A} 

In concrete NP models, the contributions to the 6 Wilson coefficients Cj \ 10 are typ- 
ically highly correlated and not all of them receive NP contributions. In addition to a 
completely generic case, we will therefore consider several restricted scenarios, that are 
each representative for a vast class of models: 

1. Real left-handed currents, d G H, C[ = 0. This is realised e.g. in models with 
MFV in the definition of [II [2j, i.e. no CP violation beyond the CKM phase. 

2. Complex left-handed currents, Cj G <D, C\ = 0. This is realised e.g. in models 
with MFV and flavour-blind phases. 

3. Complex right-handed currents, C[ G <D, Ci = 0. 

4. Generic NP, Q G <D, C[ G C. 

We remark that the results of the MFV setup analysed in |11| are recovered as a limiting 
case of our scenario 1 when the flavour mixing angles are taken to be CKM-like Q 

3.2.1. Real left-handed currents 

In this scenario, there are only 3 free parameters, the real Wilson coefficients CV^io- 
Fig.[3]shows the 68% and 95% confidence regions in the planes of two Wilson coefficient^] 
While departures from the SM point are already quite strongly restricted, a crucial 
feature is the presence of a second allowed region, characterised by a simultaneous sign 
flip of Cio as well as of the effective low-energy coefficients C| ff and C| . A sign flip of 
only Cg S and C% s or C\q and C^ s is excluded mainly by high-g 2 B — > K* ^ fjT data, 

7 In principle, in MFV, we should also account for the constraints arising fro m b — > d and s — > d 
transitions, such as T5R(Kl — > fj, + /j,~), BR(if — > ir£ + £~) or BR(B — > Xdj) [85]- In practise, the 
latter turn out to be less stringent compared to the constraints from b — > s transitions. 

8 Here and in the following, by "C.L." we mean Bayesian confidence regions, i.e. regions containing 68% 
or 95% of the Markov chain points. 
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Figure 3: Constraints between Wilson coefficients in the scenario with real left-handed 
currents. Shown are 68% and 95% C.L. regions. 



as was discussed already in 14 . Interestingly, a sign flip of only Cg S and Cio is now 
excluded as well by the precise ^4fb measurement at low q 2 (see section 3.1 ). We remark 



that an overall sign flip of all Wilson coefficients cannot be excluded by low-energy data 
alone, since all observables involve squared amplitudes or interference terms, which are 
invariant under such a sign flip. On the other hand, a NP model which generates effects 
in CV, Cg and C\q which are each twice as large as their SM contributions seems highly 
unlikely. In the region with SM-like signs, we obtain the constraints 

C 7 NP G [-0.15, 0.03] , C 9 NP G [-1.1, 1.6] , Cf p G [-1.2, 1.6] , (18) 

at 95% C.L. These constraints can be translated into bounds on an effective NP scale A 
that suppresses NP contributions to the corresponding higher dimensional operators in 
the effective Hamiltonian 

bs bs bs 

Hj = C i,O 7 + %O 9 + ^O 10 +h.c. (19) 



A 2 ' ' A2 y 1 A 



2 
10 



Assuming the coefficients c^ s to be 1, and using 95% C.L. bounds on the absolute values 
of the Wilson coefficients we obtain 

A 7 > 55 TeV , A 9 > 20 TeV , A 10 > 21 TeV . (20) 

These bounds on the effective NP scale are still weaker than the ones that can be obtained 
from considering dimension 6 operators that contribute to B s mixing |86|. 



3.2.2. Complex left-handed currents 

Complex contributions to the Wilson coefficients CV g io with vanishing C[ are predicted 
e.g. in models with MFV and flavour-blind phases. Fig. [4] shows the 68% and 95% 
confidence regions in the planes of two Wilson coefficients, omitting plots lacking a 
correlation. At 68% C.L., we observe two solutions for the real parts of the Wilson 
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Figure 4: Constraints between Wilson coefficients in the scenario with complex left 
handed currents. Shown are 68% and 95% C.L. regions. 



15 



coefficients, corresponding to a SM-like case and a case with a simultaneous sign flip 
in the low-energy values of C 10 , Cf and Cf. However, the room for NP is much 
larger than in the case without non-standard CP violation since the increased number of 
free parameters allows for compensations among different contributions. In particular, 
sizable imaginary parts are allowed for all three coefficients and this implies, in turn, 



potentially large effects in CP violating observables, as will be quantified in sec. 3.2.6 
A strong anti-correlation can be observed between the imaginary parts of Cg and Cio, 
driven mostly by the forward-backward asymmetry in B — > K*fi + fi~ at high q 2 . 

3.2.3. Complex right-handed currents 

Fig. [5] shows the 68% and 95% confidence regions in the planes of two Wilson coeffi- 
cients, omitting plots lacking a correlation. One of the most prominent differences we 
can observe compared to the case of left-handed curents, is the absence of two solutions 
for the real and imaginary parts of the Wilson coefficients. This can be mainly traced 
back to the lack of interference between right-handed currents with the SM contribu- 
tions in inclusive decays, which forbids sign flips in the low-energy values of the Wilson 
coefficients. Also in this scenario the NP room for the imaginary parts of C 7 , C' g and 
C' w is quite sizable and this will induce large effects in CP violating observables, as 



discussed in sec. 3.2.6 However, some of the correlations among low energy observables 
turn out to be different compared to the case of left-handed current, providing a tool to 
distinguish the two scenarios. 

3.2.4. Generic NP 

In this fit, we allowed all six Wilson coefficients cig 10 to assume arbitrary complex 
values. Fig. [6] shows the 68% and 95% confidence regions in the planes of two Wilson 
coefficients, omitting plots lacking a correlation. Sizable effects are not ruled out in 
any of the Wilson coefficients. However, some observations can be made, which are, by 
means of the generality of this fit, valid for any theory beyond the SM. 

• For C7, Cg and C\q there is little room left for constructive interference of real NP 
contributions with the SM. Concretely, we find at 95% C.L. 

Re C^ p > -0.01 , Re Cf p < 0.2 , Re Cft p > 0.5 . (21) 

The slight preference for non-standard Cio is again driven by the tension between 
SM and experiment in BR(f? -)• X s £ + £-) at high q 2 and F L (B -)• K*n + fj,~) at 
low q 2 . 

• a sizable negative real NP contribution to Cg requires comparably large (with 
respect to the SM) positive contributions to Cio and CV- 

• Large imaginary parts for all coefficients and large chirality-flipped coefficients are 
still allowed. 
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Figure 5: Constraints between Wilson coefficients in the scenario with complex right- 
handed currents. Shown are 68% and 95% C.L. regions. 
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BR(B u + u~) 


BR(B S t+t~) 


1 V^l/ [1,6] 1 


\(Aa)n si 1 
1 V^S/fl.B] 1 


1 ( An) n Kl I 


[1,6] 


Real LH 


[1.0,5.6] x 10~ 9 


[2, 12] x 10~ 7 


o 


o 


o 





Complex LH 


[1.0,5.4] x 10~ 9 


[2, 12] x 10~ 7 


< 0.31 


< 0.15 





o 


Complex RH 


< 5 6 x 10~ 9 


< 12 x 10~ 7 


< 0.22 


< 0.17 


< 0.12 


[-0 06 151 


Generic NP 


< 5.5 x 10~ 9 


< 12 x 10~ 7 


< 0.34 


< 0.20 


< 0.15 


[-0.11,0.18] 


LH Z peng. 


[1.4,5.5] x 10~ 9 


[3, 12] x 10- 7 


< 0.27 


< 0.14 








RH Z peng. 


< 3.8 x 10~ 9 


< 8 x 10- 7 


< 0.22 


< 0.18 


< 0.12 


[-0.03,0.18] 


Generic Z p. 


< 4.1 x 10~ 9 


< 9 x 10- 7 


< 0.28 


< 0.21 


< 0.13 


[-0.07,0.19] 


scalar current 


< 1.1 x 10~ 8 


< 1.3(2.3) x 10~ 6 















Table 3: Predictions at 95% C.L. for the branching ratios of B s — > and B s — > t + t~ 

and predictions for \ow-q 2 angular observables in B — > K*fi + fi~ (neglecting tiny 
SM effects below the percent level) in all the scenarios. The scenarios "Real 
LH" , "Complex LH" , "Complex RH" , "Generic NP" , "LH Z peng." , "RH Z 



peng.", and "Generic Z p." correspond to the scenarios discussed in sec. |3.2.1 
3.2.3| sec. |3.2.4[ sec. |4.1.1[ sec. |4.1.2[ and sec. |4.1.3[ respectively, 



sec. 



3.2.2 



sec. 



assuming negligible (pseudo) scalar currents. In the scenario "scalar current" 
only scalar currents are considered. The number quoted for B s —> t + t~ in the 
"scalar current" scenario refers to the maximum value for its branching ratio 
in the case of dominant scalar (pseudoscalar) currents. 



The last point highlights the importance of measuring observables sensitive to right- 
handed currents and to CP violation, such as the B —> K* fi + /-i~ observables ^7,8,9 and 
S3. 

3.2.5. Predictions for B s — >■ /i -1- /! - and B s — > t + t~ 

Due to its sensitivity to scalar currents, we did not include the branching ratios of 
B s — > fi + ^L~ and B s — > t + t~ in our fits. Instead, using the constraints on cj'J, we can 
now give upper limits on the branching ratios in the considered scenarios assuming the 
scalar and pseudoscalar Wilson coefficients to be negligible. These bounds are useful 
since an observed violation of them would imply the presence of scalar currents. The 
values we find at 95% C.L. are listed in table [3j In the scenario with right-handed 
currents as well as in the case of generic NP, a significant suppression of the branching 
ratios is possible. The upper bounds on BR(i? s — > in all cases are approximately 

a factor of 2 below the current experimental bound Q and correspond roughly to a 50% 
enhancement of the branching ratio with respect to the SM. Due to the new measurement 
of B — > K*fi + [i~ angular observables, they are stronger than similar bounds presented 



in the literature before 14 . 
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We stress that scalar current effects in B s — > /^ + ^ _ could still enhance the branching 
ratio over its current experimental bound. 

Let us also mention that in the case of B d — > both the effects induced by 

the corresponding b — )■ d Wilson coefficient c['q and scalar currents can enhance the 
branching ratio of B d — > /U + /U _ over the experimental bound since the corresponding 
constraints from b — > d£ + £~ processes are weaker. In addition, we remark that the ratio 
BR(i? c / —7- fi + /x~)/BR(i? s —7- can significantly depart from the SM as well as the 

MFV predictions BR(B d -> n + n~)/BR(B s « \ V td /V t s\ 2 in both directions. 

Finally, we discuss the allowed values for the branching ratio of B s — > t + t~ . From 
the general expressions of BR(I3 S — > in presence of NP (see eq. [5]) one has 

BR(fl.->rV) ^ / Aml\ l ' 2 ml (l - 4m*/m%.) \S\ 2 + \P\ 2 

BR(£ S ~ I m| J m 2 \S\ 2 + \P\ 2 ' 1 j 



where S and P have been defined in eq. [6j 

les the dom: 

BR{B S r+T 



In the case where cj'J provides the dominant NP effects, one obtains 



212 , (23) 



BR{B S -»■ /z+/z~ 

which implies, in particular, the SM prediction for the branching ratio of B s — > t + t~ 

BR{B S -»■ t+Osm = (7.7 ± 0.8) x 10~ 7 . (24) 

In the case where (pseudo)scalar current effects dominate, BR(£? S — > /i + /i~) can saturate 
the current experimental bound while for BR(i? s — > t + t~) we get 

. BR(B S -> r+T") ^ , , 

120< R ' + _ <212, 25 
BR(£> S — )• //^/i ) 

where the lower (upper) bound in eq. [25] correspond to the case where the scalar (pseu- 
doscalar) contribution dominates. 



Combining eq. (|3j) with eq. ( |25| , it turns out that BR(f? s — > t + t~) < 2 x 10~ b (see 
also table [3]) and it has to be seen whether such values might be within the reach of 
LHCb. However, we stress that this upper bound relies on the assumption that the 
(pseudo)scalar Wilson coefficients C^ p are linearly proportional to the lepton Yukawa 
couplings, as discussed in sec. [2] If we relax this assumption, as it might be the case in 
models like R-parity violating SUSY |87| or models with enhanced couplings with the 
third lepton generation (compared to the linear scaling assumed throughout this work), 
BR(i? s — > t + t~) could get in principle much larger values than 10 -6 . 

3.2.6. Predictions for B K*fj, + V~ 

Figure [7] shows the predictions for the T-odd B — > K*[i + fi~ CP asymmetries A7 and 
A$ at low q 2 for the scenarios with complex left-handed currents, complex right-handed 
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Figure 7: Fit predictions for the low-g 2 CP asymmetries (^7,8) m B — > K*fi + fi~ in 
the case of complex left-handed currents (left), complex right-handed currents 
(centre) and generic NP (right). Shown are 68% and 95% C.L. regions. 



currents and for generic NP. In the absence of right-handed currents, one finds an anti- 
correlation between Aj and ^4s which has already been found in models where only C7 
contributes (55 , 82 , 84 (see also [88]), but is shown here to hold under more general con- 
ditions. At 68% C.L., one finds a preference for non-standard CP asymmetries driven 
mostly by the tension between SM and experiment in Fl(B — > K* n + n~) at low q 2 . Sim- 
ilarly, in the absence of complex left-handed currents, one finds an opposite correlation. 
In the generic case, there is no correlation at all. Interestingly, in all three scenarios, 
large effects in both asymmetries are still allowed, with the numerical bounds listed in 
table [3| Future measurements of Aj and A$ at LHCb will thus be crucial to constrain 
the imaginary parts of the Wilson coefficients entering the B — > K* ji + fi~ decay. 

Also shown in table [3] are the predictions for the CP asymmetry Ag and the CP- 
averaged angular coefficient S3 at low q 2 , both of which are tiny in the SM but can be 
sizable in presence of right-handed currents. Indeed, both observables can assume values 
in excess of 10% in the complex right-handed scenario and for generic NP. 

We note that the CDF measurement of S3 and Ag shown in the last two rows of table [3] 
currently puts no significant constraints on NP, yet. Future measurements at LHCb with 
errors of the order of 0.1 will however put important constraints on CP-violating or CP- 
conserving right-handed currents. 



4. Analysis of flavour-changing Z couplings 



Tree level FCNC couplings of the Z boson can appear in a number of NP scenar- 
ios. Prominent examples are the SM with four non-sequential generations of quarks 
and models with an extra U(l) symmetry 89,90 . Moreover, the contributions to the 



semi-leptonic operators are dominated by Z penguins, i.e. loop-induced modified Z cou- 
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plinga I in many theories, e.g. in the MSSM (92pi| . It is therefore interesting to consider 
the effects in a framework with modified sbZ couplings, which can be parametrised by 



the effective Lagrangian 92 



£%f = - -^-^m 2 z c w s w V; b V ts Z» (Z L s 7fl P L b + Z R s^P R b) , (26) 
with s,„ = sin6>,„ and c,„ = cos 6,„. In this class of models one finds 



c io - Zf v , C' w - Z R , (27) 

C7 9 NP = -Zf(l - 44) , C 9 = -Z R (1 - 44) . (28) 

The contributions to Cg are strongly suppressed by the small vector coupling of the Z 
to charged leptons (1 — 44) ~ 0.08. 

The modified sbZ couplings also modify decays with a neutrino pair in the final state, 
so one obtains correlations between b — > svv and b —> s£ + £~~ observables. Writing the 
b —7- svv effective Hamiltonian as 

n cS = --j£-V tb V t * s {C l O l + C r O r ) + h.c. , (29) 

e 2 

Ol,r = ^ i^ P L,Rb) {v^Plv) , (30) 



The effective Lagrangian (26) leads to 

C L = Cf i + Zf ) Cr = Z r , (31) 



where Cf M = -6.38 ± 0.06 [93 

Finally, an effective tree level contribution to B s -B s mixing is generated by the ex- 
change of a Z with modified sbZ coupling. Its contribution to the mixing amplitude can 
be written at the scale mz as [93] 



{B s \n\B s y bZ _ 4o4 2 2 

(B.\H\B.)™ ~ ^SoK) (Zl ~ 3 ' 5 ZLZR + Zr) ■ (32) 
We consider three scenarios in the following: 

• Left-handed modified Z couplings, Z L G <D, Cf p G C, Z R = C' 7 = 0, 

• Right-handed modified Z couplings, Z R G C, C 7 G C, Z^ p = p = 0, 

• Generic modified Z couplings, Zi >R G C, C^p G C, 

allowing for non-standard CP violation in all cases. The generic case covers all NP models 
where contributions to the semileptonic operators are dominated by Z penguins. This 
includes in particular the general MSSM (92 , 94 . The fitting procedure is as described 
in section [3T2l 
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Figure 8: Constraints on C7 and the modified Z coupling in the scenario with left-handed 
couplings only. Shown are 68% and 95% C.L. regions. 



4.1. Constraints on modified Z couplings 
4.1.1. Left-handed modified Z couplings 

Fig. [8] shows the 68% and 95% confidence regions in the complex planes of C7 and Z^ 
as well as the correlation between the real parts of C7 and Z^. The constraint in the 
C7 plane is basically identical to the constraint in the absence of semileptonic operators 
shown in the upper left plot of figure [2j which is in contrast to the corresponding con- 
straint in the presence of Cg and C\q shown in figure [4j where large effects in Re(C7) 
(and in particular a sign flip) were allowed. This can be traced back to the suppression 
of Cg compared to C10 in the Z penguin scenario, making impossible a simultaneous 
sign flip of Cj S and C| at low energies, which would be required in particular to meet 
the constraint from Af&(B — > K*{j, + [i~) at low q 2 . For the same reason, the correlation 
of the real parts of C7 and Zl only shows one solution. 



4.1.2. Right-handed modified Z couplings 

Fig. [9] shows the 68% and 95% confidence regions in the complex planes of C 7 and Zr 
as well as the correlation between the real parts of C' 7 and Zr. The constraint in the C' 7 
plane is very similar to the constraint in the absence of semileptonic operators shown in 
the corresponding plot of figure [2j The negative values preferred for the real part of the 
right-handed Z coupling, i.e. for Re(C( ), arises from low- and high q 2 B — >■ K*[i + [i~ 
data, as can be seen in the corresponding plot of fig. [2} 



4.1.3. Generic modified Z couplings 



Fig. 10 shows the 68% and 95% confidence regions in the complex planes of the Wilson 
coefficients C 7 and the couplings Zl : r in the case of generic modified Z penguins. While 



9 We remark that the distinction between Z penguins and other contributions is in general gauge 
dependent. However, this gauge dependence is weak if Z penguins dominate 



91 



9:3 
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Figure 9: Constraints on C-j and the modified Z coupling in the scenario with right- 
handed couplings only. Shown are 68% and 95% C.L. regions. 
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Figure 10: Constraints on Cn and the modified Z couplings in the scenario with generic 
modified Z couplings. Shown are 68% and 95% C.L. regions. 
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Figure 11: Fit predictions for the low-q 2 CP asymmetries (^4.7,s) m B — > K*fi + fi~ for 
the scenario with left-handed (left), right-handed (centre) or generic (right) 
modified Z couplings. Shown are 68% and 95% C.L. regions. 



the room for NP is larger than in the more constrained previous cases, also in the generic 
there are no disjoint solutions for the Wilson coefficients. We are thus lead to conclude 
on a model-independent basis that if the NP contributions to semi-leptonic operators 
are dominated by Z penguins, the real parts of the Wilson coefficients 67,9,10 at low 
energies must have the same sign as in the SM. 



4.2. Fit predictions 



4.2.1. Predictions for B„ 



T + T - 



and B 



K' 



Analogously to section 3.2.5, we can give fit predictions for BR(5 S — > and 
BR(i? s — > t + t~) in the absence of scalar currents and for B — > K*/i + fi~ observables in 
the considered modified Z coupling scenarios based on the constraints obtained in the 
global fit. The allowed ranges are shown in table [3} In the generic case, the preference 
for smaller values of the B s — > li \i~ and B s — > t + t~ branching ratios is due to the 



negative values preferred for Re(Zfj) (cf. section [4.1. 2[ ), i.e. for Re(Cj ), which leads to 
a destructive interference with the SM in the decay amplitudes, see eq. Q. 



Figure 11 shows the prediction for the B — > K*fi + fi CP asymmetries A-j and A% 
at low q 2 in all three scenarios. The predictions are similar to the corresponding ones 
obtained for generic Cg\ shown in figure 7| so the comments made there apply here as 
well. Also for the observables S3 and Ag we find predictions that are similar to the cases 
discussed in section 13.2.61 The results are summarised in table [3l 



4.2.2. Predictions for B s mixing 

Since the real and imaginary parts of the left- and right-handed Z couplings are con- 
strained by b — > s£ + £~ processes not to be significantly larger than the SM value of 
the (real) left-handed Z coupling, the Z exchange contribution to B s mixing, which is 
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negligible in the SM, cannot lead to sizable deviations from the SM. Concretely, in the 
considered scenarios we find, at 95% C.L., 



(33) 
(34) 
(35) 

Such NP contributions are well within the range allowed by the measurement of the B s 



left-handed mod. Z couplings: 




qSMi 

D ij)(f> \ 


< 


0.008 


right-handed mod. Z couplings: 




cSMi 


< 


0.014 


generic mod. Z couplings: 




cSMi 
D ijxf> \ 


< 


0.04. 
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4.2.3. Predictions for b suv decays 

The two exclusive b — > svv decays, B — > (K,K*)vi>, and the inclusive one B -4 X s vv 
give access to four observables sensitive to NP: the three branching ratios and the K* 
longitudinal polarisation fraction Fl in B — > K*vv. However, the observables are not 
all independent since they depend on only two real combinations of the complex Wilson 



coefficients Cl and Cr 93 96 97 



^\C L \ 2 + \C R \ 2 _ -Re(C L C* R ) 



KcysM| ' ' \c L \* + \c R \*- 

For the central values of the hadronic parameters, one obtain 



(36) 



10 
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BR(B->K*vD) = 6.8 x 10~ 6 (1 + 1.31 r])e 2 , (37) 
BR(B -4 Kuv) = 4.5 x 10" 6 (1-2 r?)e 2 , (38) 
BR(B -»• X a vv) = 2.7 x 10" 5 (1 + 0.09 r])e 2 , (39) 

{Fl){B — )■ K*vv) = 0-54 ^If^ , (40) 

where {Fl) refers to the ratio of the branching ratio into a longitudinal K* over the total 
branching ratio. It can be extracted from the angular distribution of the K* — > Kir decay 
products. 

In the scenario with left-handed modified Z couplings only, one has r\ = 0, Fl is 
SM-like and all the branching ratios are merely scaled by a common factor. We obtain, 
at 95% C.L., 

g2 = BR(B -> K*vv) = BR(B -> Kvv) 

RR(B -4 K*vv)su BR(B -4 Kvv) S u ' ' ' ' 

In the scenario with right-handed modified Z couplings and SM-like left-handed cou- 
plings, one finds an anticorrelation between the two experimentally most promising 

10 A lower central value for BR(B — > Kvv) is obtained if the experimental value of BR(_B — > K£ + £~) is 
used, assuming the latter decay to be SM-like [98] . Here we allow both decays to deviate from the 
SM prediction. We treat the B — > r(— > Kv)v contribution [99] as a background to be subtracted 
from the experimental result. 
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2 4 6 8 10 12 14 



10 6 xBR(B^£vv) 

Figure 12: Fit prediction for the branching ratios of B — > K^*'vv for generic modified Z 
couplings. Shown are 68% and 95% C.L. regions. 



modes, B -»• Kvv and B -)■ K*vv (see ppjoj ) . At 95% C.L., we find 



BR(B -> K*vv) . , BR(B -»• #i/z?) . . (F L ) . , 

^ ^- G 0.6, 1.0 , H v^- S 1.0, 2.6 , / y G 0.4, 1.0 , 

BR{B^K*vv) SM 1 ' J ' BR(i?^K^)sM L ' J ' (F L ) SU 1 ' J ' 

(42) 

The preference for a suppression of Fl and BR(i3 — > K*vv) but an enhancement of 
BR(i? —7- Kvv) is again due to the negative values preferred for Zr commented on in 
section 4.1.2 The large enhancement possible for BR(£? —> Kvv) is close to the current 
experimental bound BR(B + -> K+vv) < 13 x 10" 6 flOlj . 

In the case of left- and right-handed Z couplings, the correlation between the decays 
carries information on the size of left- vs. right-handed currents and is thus a valuable 



probe of the chirality structure of NP. In figure 12, we show the fit prediction for B — > 
Kvv and B — > K*vv. We observe that the SM point is allowed at about 68% C.L. 
Also in this case, only a small enhancement or a sizable suppression is allowed for the 
B —7- K*vv branching ratio, while an enhancement of the B — > Kvv branching ratio up 
to a factor of 3 is possible. For Fl, we find at 95% C.L. 

^ G [0.4, 1.1] . (43) 



(Fl)su 

In all the scenarios, the full allowed range of branching ratios can be probed at the 
next-generation B factories (38) . 

Finally, we stress that these conclusions are only valid under the assumptions that 
modified Z couplings dominate the NP contributions to b — > svv and b —> s£ + £~ semi- 
leptonic operators. Much larger effects in B — > K^vv are possible in models where this 
is not the case, e.g. models with a non-universal Z' coupling stronger to neutrinos than 
to charged leptons |93|. 
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5. Conclusions 



Rare decays with a b — > s transition offer excellent opportunities to probe the flavour 
sectors of extensions of the SM. The effects of new heavy degrees of freedom in these 
processes can be parametrised by modifications of the Wilson coefficients of local, non- 
renormalizable operators, which allows to constrain such NP effects in a model-independent 
way. In this work, we analysed the constraints on the Wilson coefficients that follow from 
the currently available experimental data on b — > s rare decays. We took into account 
the measurements from the B factories of the branching ratio of the radiative B — > X s j 
decay, of the time-dependent CP asymmetry in B — > K*^, of the branching ratio of the 
inclusive B — > X s £ + £~ decay, Belle and CDF data on the branching ratio and angular 
distribution of the exclusive B —> K* ji + decay and in particular the recent LHCb 
results on the branching ratio and angular distribution of B — > K* : // + fj,~ . 

The constraints on the Wilson coefficients are obtained by using a \ 2 function, which 
depends on the Wilson coefficients and contains the theory predictions for the observables 
and experimental averages as well as the corresponding uncertainties. 

We have analysed the following scenarios where: 

1. the dominant NP effects are captured already by one complex Wilson coefficient 
or by a pair of real Wilson coefficients. This is a representative case of many 



NP models like the MSSM with MFV and flavour blind phases 81,82 , non-MFV 



SUSY models 83,84, 1021 and also models with dominance of Z penguins; 



2. the NP effects are accounted for by means of the full set of the 6 complex Wilson 
coefficients C 7 'g 10 . 

3. the dominant NP effects in the semi-leptonic operators arise from non-standard 
flavour changing Z couplings. 

While we refer to sections [3] and [4] for a detailed description of all our results, we want 
to emphasise here the following main messages: 

• At the 95% C.L., all best fit regions are compatible with the SM. 



The combination of inclusive and exclusive b — > sv t~ observables exclude sign 
flips in various low-energy WCs. That is, the SM is likely to provide the dominant 
effects in low energy observables. In particular, we show that 

— sign flips in CV, Cg or C\q are excluded if NP enters dominantly through Z 
penguins, 

— only a simultaneous sign flip of C7, C9 and C10, which cannot be excluded by 
low-energy data alone, is allowed in the absence of non-standard CP violation 
or right-handed currents. 

The new A-p-Q measurement at low q 2 from LHCb is already quite effective in 
constraining NP effects. The same is true for the time-dependent CP asymmetry 
in B K*j. 
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• High-g 2 data onB-> K*[i + fi are competitive with and complementary to the 
low-q 2 ones. 

Moreover, we have investigated the implications of the above constraints and the future 
prospects for observables in b — > sl + l~ and b — > svv transitions in view of improved 
measurements. In particular, we find that 

• in the presence of non-standard CP violation, the low-g 2 angular CP asymmetries 
Ai and As in B — > K*fj, + fi~ can reach up to ±35% and ±20%, respectively. 

• in the presence of right-handed currents, the low-g 2 angular observables Ag and S3 
in B — > K*fi + /j,~ can reach up to ±15%. 

• in the absence of (pseudo) scalar currents, BR(£? S — > and BR(1? S — > t + t~) 
can be enhanced at most by 50% over their SM values, mainly due to the new 
measurement of B — > K*fi + fi~ angular observables. In contrast, if scalar currents 
are at work, BR(B S — > can still saturate the current experimental bound 
and BR(1? S — > t + t~) can be enhanced by a factor of 3. 

• if NP in b —> s£ + £~ and b — > svv processes is dominated by left-handed Z penguins, 
an enhancement of the branching ratios of B — > K^vv by more than 30% is 
unlikely. If right-handed Z penguins are present, B — > Kvv can saturate the 
present experimental bound, while B — > K*vv is unlikely to be enhanced. 

The first two points highlight the importance of measuring observables in the B — > 
K* n + /i~ angular distribution sensitive to right-handed currents and CP violation. Such 
measurements would be crucial to lift degeneracies in the space of Wilson coefficients 
which make it difficult at present to put strong constraints on individual coefficients in 
a completely generic NP model, as our analysis in sections 3.2.2 and 3.2.4| showed. 



In conclusion, the present work updates and generalises previous studies providing, at 
the same time, a useful tool to test the flavour structure of any theory beyond the SM. 
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A. Statistical method 

Here we give some details on our statistical method used to obtain the constraints on the 
Wilson coefficients in section 3.2 and|4| We use a Bayesian approach with the likelihood 
function 

L(C) = e -* 2 (^)/2 , (44) 
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where C is a 12-dimensional vector containing the real and imaginary parts of the 6 
Wilson coefficients 10 and the x 2 function has been defined at the beginning of 
section [3j We sample the posterior probability distribution, defined according to Bayes' 
theorem, 

P(C)= (45) 
j L{C')it{C')dO 

by means of a Markov Chain Monte Carlo (MCMC) analysis using the Metropolis- 
Hastings algorithm (see e.g. the review on statistics in 103] ). As proposal density, 



we use a multivariate Gaussian, whose width is optimised to tune the acceptance rate. 
We use a flat prior, vr(C) = 1, in the general case and a multivariate 5 function in the 
restricted scenarios. The stationary density distribution of points in the Markov chain 
is proportional to P(C). Constraints on two-dimensional subspaces are obtained by 
projecting the points onto this plane. Predictions for observables presented in sections |3.2| 
and [4] are obtained by evaluating the observable for central values of the theoretical 
input parameters at each point in the chain and interpreting their density distribution 
as posterior probability for the observable. 

Two-dimensional confidence regions in sections 3.2 and [4] are obtained by determining 
contours of constant posterior probability which contain 68% (or 95%) of the Markov 
chain points. Analogously, one-dimensional confidence regions are highest posterior den- 
sity intervals, i.e. the posterior probability is higher everywhere inside the interval than 
outside, and they contain 68% (or 95%) of the points. 



B. Effective Wilson coefficients 

In sections [3] and [4j we have put constraints on NP contributions to the Wilson co- 
efficients at a matching scale of 160 GeV. Here we give the relation to the effective 
Wilson coefficients at low energies, which are the quantities relevant for the evaluation 
of observables. 

In low-energy observables, the coefficients Cj and Cg always appear in a particular 
combination with four-quark operators (which can be found e.g. in [l9]) in matrix 
elements. It hence proves convenient to define effective coefficients Cf^, which are given 



by 1 104 



Cf = C 7 -±C 3 -^C 4 -fc 5 -yC 6) (46) 
Cf = C 9 + Y(q 2 ), (47) 
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with 

Y(q 2 ) = h(q 2 , m c ) f | d + C 2 + 6C3 + 60C 5 

- \ Hq 2 , m b ) (jC 3 + | C 4 + 76C 5 + y C 6 

1 o / 4 64 \ 4 64 64 

--% 2 ,0) fC3 + -C 4 + 16C 5 + -C 6 J +-C 3 + -C 5 + -C7 6 (48) 

at leading order in a s and Aqcd/^-6- Beyond the leading order, there are perturbative 
corrections as well as power corrections, which differ for inclusive and exclusive b —> 
s£ + £~ decays. We refer the reader to Refs. [44|[45|[47| - [5T|[69|pf4] for these corrections. 
Ref. pTO] contains the expressions for the doubly Cabibbo-suppressed contribution to 



(46)— (47) which is relevant for the SM prediction of the B — > K* fJ, + fJ, CP asymmetries. 



In the SM, one finds at the scale \x b = 4.8 GeV, to NNLL accuracy |55| , 

Cf(M&) = -0.304, CM =4.211, do(M6) = -4.103, (49) 

while the primed coefficients are negligible. Beyond the SM, but assuming the four-quark 
operators to be free of NP, one has 

Cf(p b ) = Cf ' SM ( W ) + C 7 NP (^) , (50) 

Of (14) = Cf > SU (» b ) + C 9 NP , (51) 

Cio = C 1Q + Cf p , (52) 

C> 7 ( f x b ) = C™ p ( f x b ), (53) 

C' 9 ,M = C'* P . (54) 

While the NP contributions to Cq} and c[q do not run, C^p do and they mix with C^ 
under renormalization. With leading order running, as is appropriate for NP contribu- 
tions evaluated at one loop, from a high matching scale fih = 160 GeV, one finds 

cf NP (/i fe ) = 0.623 cf )NP (/i h ) + 0.101 cP NP (fi h ) . (55) 

Since the low-energy observables are sensitive to C) (/•*&), the constraints we presented 
in sections [3] and [4] on C^ ^ NP (/i^) for vanishing Cg ^ NP can be interpreted as constraints 
on Cf )NP Gu h ) + 0.162 cP NP (nh) for non-standard cf } . 
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